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ABSTRACT 

We first present the results of numerical simulations on formation processes and 
physical properties of old globular clusters (GCs) located within clusters of galaxies 
("intracluster GCs") and in between clusters of galaxies ( "intercluster GCs"). Our 
high-resolution cosmological simulations with models of GC formation at high red- 
shifts (z > 6) show that about 30 % of all GCs in a rich cluster can be ragarded as 
intracluster GCs that can freely drift being trapped by gravitational potential of the 
cluster rather than by the cluster member galaxies. The radial surface density profiles 
of the simulated intracluster GCs are highly likely to be flatter than those of GCs 
within cluster member galaxies. We also find that about 1% of all GCs formed before 
z > 6 are not located within any virialized halos and can be regarded as "intercluster" 
(or "intergalactic" ) GCs. We discuss the dependences of physical properties of intr- 
acluster and intercluster GCs on the initial density profiles of GCs within low-mass 
dark matter halos at high redshifts (z > 6). 
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1 INTRODUCTION 

Recent observational studies of globular clusters (GCs) in 
clusters of galaxies have suggested that there can be a pop- 
ulation of GCs that are bounded by cluster gravitational 
potentials rather than those of cluster member galaxies 
(e.g., West et al. 1995; Bassino et al. 2002, 2003; Jordan 
et al. 2003), though the existence of these intracluster 
GCs (ICGCs) in the Coma cluster is observationally sug- 
gested to be highly unlikely (e.g., Marin- Franch & Apari- 
cio 2003). Structural and kinematical studies of a popula- 
tion of very bright star clusters — known as "ultra-compact 
dwarfs (UCDs)" — have also suggested that these clusters 
can be also freely floating intracluster objects (Mieske et al. 
2004; Drinkwater et al. 2005; Jones et al. 2005). Physical 
properties of intracluster stellar objects such as ICGCs and 
PNe are considered to be sensitive to dark matter properties 
and cluster- related physical processes (e.g., tidal stripping of 
GCs and hierarchical growth of clusters) and thus provide 
some fossil information on formation of galaxies and clusters 
of galaxies (e.g., Arnaboldi 2004 for a recent review). 

Although there have been developments on the observa- 
tional front, there has been little theoretical and numerical 
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works carried out as to how ICGCs are formed in clusters en- 
vironments (e.g., Forte et al. 1982; Muzzio et al. 1987; Bekki 
et al. 2003) . These previous models showed that tidal strip- 
ping of GCs from cluster member galaxies though galaxy- 
galaxy and galaxy-cluster interaction is a mechanism for 
ICGC formation. These previous works however used fixed 
gravitational potentials of already virialized clusters and 
accordingly could not discuss how ICGCs in a cluster are 
formed as the cluster grows through hierarchical merging of 
smaller groups and clusters. Thus it remains unclear (1) how 
ICGC are formed under the currently favored cold dark mat- 
ter (CDM) theory of galaxy formation and (2) what physical 
properties ICGCs can have if their formation is closely as- 
sociated with hierarchical formation of clusters. 

The purpose of this Letter is thus to demonstrate, for 
the first time, how ICGCs are formed during hierarchical 
formation of clusters of galaxies, based on high-resolution 
cosmological simulations that can follow both hierarchical 
growth of clusters through merging of smaller subhalos and 
dynamical evolution of old GCs. We also discuss physical 
properties of GCs that were formed within subhalos at z > 6 
yet are not within any virialized halos at z — 0: These GCs 
can be regarded as "intergalactic" (van den Bergh 1958) or 
"intercluster" GCs. For convenience and clarity, GCs within 
any virialized halos with the masses larger than 3 x 1O 9 M0 
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Figure 1. Left: The large scale distributions of INGCs (magenta big dots) and HGCs (cyan dots) projected onto the x-y plane at z = 
in the fiducial model. Here INGCs represent GCs that are not within any virialized halos at z = and thus include intergalactic and 
intercluster GCs. HGCs are GCs that are within virialized halos at z = 0. Right: Distributions of GCs projected onto the x-y plane at 
2 = for a cluster-scale halo with the total mass of 6.5 X 1O 14 M0. GCs within circles represent those within tidal radii of galaxy-scale 
halos and the radii of the circles represent the tidal radii. GCs that are not within any circles are regarded as ICGCs. 



are refereed to as Halo Globular Clusters (HGCs) and the 
meanings of these acronym are given in Table 1. 



2 THE MODEL 

We simulate the large scale structure of GCs in a ACDM 
Universe with Q = 0.3, A = 0.7, H = 70 km s" 1 Mpc" 1 , 
and as = 0.9 by using the Adaptive Mesh Refinement 
N— body code developed by Yahagi (2005) and Yahagi et al. 
(2004), which is a vectorized and parallelized version of the 
code described in Yahagi & Yoshii (2001). We use 512 3 colli- 
sionless dark matter (DM) particles in a simulation with the 
box size of 70/i _1 Mpc and the total mass of 4.08 x 10 16 M Q . 
We start simulations at z = 41 and follow it till z = in 
order to investigate physical properties of old GCs outside 
and inside virialized dark matter halos at z = 0. We used the 
COSMICS (Cosmological Initial Conditions and Microwave 
Anisotropy Codes), which is a package of fortran programs 
for generating Gaussian random initial conditions for non- 
linear structure formation simulations (Bertschinger 1995, 
2001). 

The way of investigating GC properties is described as 
follows. Firstly we select virialized dark matter subhalos at 
z — Zform by using the friends-of-friends (FoF) algorithm 
(Davis et al. 1985) with a fixed linking length of 0.2 times 
the mean DM particle separation. The minimum particle 
number N m i n for halos is set to be 10. For each individ- 
ual virialized subhalo with the half-mass radius of Rh , some 
fraction (/ gc ) of particles within -Rh/3 are labeled as "GC" 
particles. This procedure for defining GC particles is based 
on the assumption that energy dissipation via radiative cool- 



ing allows baryon to fall into the deepest potential well of 
dark-matter halos and finally to be converted into GCs. The 
value of the truncation radius (-Rtr, gc = -Rh/3) is chosen, 
because the size of the very old GCs in the Galactic GC 
system (i.e., the radius within which most Galactic old GCs 
are located) is similar to -Rh/3 of the dark matter halo in 
the dynamical model of the Galaxy (Bekki et al. 2005). We 
assume that old, metal-poor globular cluster formation is 
truncated after z = Zform, because previous theoretical stud- 
ies demonstrated that such truncation of GC formation by 
some physical mechanisms (e.g., reionization) is necessary 
for explaining the color bimodality of GCs, very high spe- 
cific frequency (<Sn) in cluster Es, and structural properties 
of the Galactic old stars and GCs (e.g., Beasley et al. 2002; 
Santos 2003; Bekki 2005; Bekki & Chiba 2005). 

Secondly we follow the dynamical evolution of GC par- 
ticles formed before z = Zf orm till z = and thereby de- 
rive locations {(x,y,z)) and velocities ((« x , v y , v z )) of GCs 
at z = 0. We then identify virialized halos at z = with the 
FoF algorithm and investigate whether each of GCs is within 
the halos. If a GC is found to be within a halo, the mass 
of the host halo (Mh) and the distance of the GC from the 
center of the halo (-R gc ) are investigated. If a GC is not in 
any halos, it is regarded as an INGC and the distance (-R nc i) 
between the INGC and the nearest neighbor halo and the 
mass of the halo (Mh,nei) are investigated. If a GC is found to 
be within a cluster-size halo (Mh > 1O 14 M0), we investigate 
which galaxy-scale halo in the cluster-scale halo contains the 
GC. We examine local mass densities around particles in a 
cluster and thereby select galaxy-scale halos that have high 
densities enough to be identified as galaxy-scale halos (BY). 
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Table 1. Meaning of acronym 

INGC INtcrgalactic GCs 

ICGC IntraCluster GCs 

GGC Galactic GCs 

HGC Halo GCs 



If we find the GC within the tidal radius (R t ) of one of 
galaxy-scale halos in the cluster-scale halo, it is regarded as 
a galactic GC (GGC): Otherwise it is regarded as an ICGC. 
-Rt is assumed to be the radius where the slope a in the GCS 
density profile of p(r) oc r a in a galaxy-scale halo becomes 
smaller than -5 (i.e., much steeper than the outer profile of 
the dark matter halo). 

Thus, the present simulations enable us to investigate 
physical properties only for old GCS owing to the adopted 
assumptions of collisionless simulations: Physical properties 
of metal-rich GCs lately formed during secondary dissipative 
galaxy merger events at lower redshifts (e.g., Ashman & Zepf 
1992) can not be predicted by this study. We present the 
results of the model with z Iorm = 6, and the dependences of 
the results on Zf OIln will be given in Bekki & Yahagi (2005, 
BY). If zf orm is closely associated with the completion of 
cosmic reionization, Zf orm may well range from 6 (Fan et 
al. 2003) to 20 (Kogut et al. 2003). Physical properties of 
hypothetical GC particles for ICGCs in a rich cluster with 
M h = 6.5 x 1O 14 M are described for the model with / gc = 
0.2 in which the number ratio of GC particles to all particles 
is 1.5 x 10~ 3 at z — 0. We adopt / gc = 0.2 so that typical 
subhalos at z = 6 can contain at least one GC particle. The 
present results does not depend on / gc at all as long as / gc > 
0.1. Physical properties of ICGCs in groups and clusters 
with different masses will be given in our forthcoming papers 
(BY). 

We assume that the initial radial (r) profiles of GCSs 
(p(r)) in subhalos at z = 6 are the same as those of the simu- 
lated dark matter halos that can have the universal "NFW" 
profiles (Navarro, Frenk, & White 1996) with p(r) oc r~ 3 in 
their outer parts. The mean mass of subhalos at z=6 in the 
present simulations is roughly 1.8 x 1O 1O M0, which is simi- 
lar to the total mass of bright dwarf galaxies. Minniti et al. 
(1996) found that the projected (R) density profiles of GCSs 
in dwarfs is approximated as p(R) oc R~ 2 , which is trans- 
lated roughly as p(r) cc r~ 3 by using a canonical conver- 
sion formula from p(R) into p(r) (Binney & Tremaine 1987). 
Therefore, the above assumption on p(r) can be regarded as 
reasonable. Thus we mainly show the fiducial model with 
p(r) similar to the NFW profiles and -Rtr, gc = Rh/3. 

Although we base our GC models on observational re- 
sults of GCSs at z—0, we can not confirm whether the above 
p(r) and i?tr, gc of the fiducial model are really the most prob- 
able (and the best) for GCSs for low-mass subhalos at z—6 
owing to the lack of observational studies of GCSs at high 
redshifts. Therefore we investigate how the numerical results 
depend on initial p(r) and -Rtr, gc of subhalos at z = 6. Since 
the dependences on p(r) and -Rtr, g c are given in details by 
BY, we briefly describe the dependence on Rti, gc , which is 
the most important dependence for physical properties of 
GCSs at z = (BY). 
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Figure 2. Projected radial density profiles (Eqc(R)) of ICGCs 
(thick) and GGCs (thin). Here GGCs represent "galactic GCs" 
that are GCs within any galaxy-scale halos in a cluster-scale 
halo. Accordingly GGCs are not GCs of the Milky Way: We here 
use this term in order to distinguish these GC populations (i.e., 
GGCs) from ICGCs. For clarity, Egc (R) normalized to the max- 
imum value in each GC population is shown. For comparison, 
the power-law density profiles with T,qq(R) oc R a for a = —1.0, 

— 1.5, —2.0, and —2.5 are shown by dotted lines. 

3 RESULT 
3.1 ICGCs 

Figure 1 shows the large scale (~ 100 Mpc) structure of 
INGCs and HGCs and the distributions of GCs in a halo 
with the total mass (Mh) of 6.5 x 1O 14 M0 corresponding to 
a rich cluster of galaxies at z = 0. It is clear from Figure 1 
that there exists ICGCs that are not bounded in any cluster 
member galaxy-scale halos, though most GCs are within the 
galaxy-scale halos. About 29 % of all GCs in this cluster can 
be classified as ICGCs with the number fraction of ICGCs 
ranging from 0.28 at R c \ < 1 Mpc (where R c \ is the distance 
between a GC and the center of the cluster) to 0.35 at 1 < 
R c i < 2 Mpc. Although the number fraction of ICGCs in the 
central 200 kpc of the cluster is only 0.02, the presence of 
such central ICGCs may well support the scenario by West 
et al. (1995) that very high 5n of ICGCs in the central giant 
Es in some clusters can be due to ICGCs. 

Figure 2 shows that the projected number density distri- 
butions (Egc) within the central few hundreds kpc are quite 
different between ICGCs and GCs within any galaxy-scale 
halos in the cluster ("galactic GCs" referred to as "GGCs") 
in the sense that Egc is significantly flatter in ICGCs than in 
GGCs. This is because most GCs in the central few hundreds 
kpc of the cluster can be identified as GCs within galaxy- 
scale halos (i.e., smaller number of ICGCs). If GCs that are 
freely drifting under the influence of the cluster potential 
are located close to the central giant halo(s) of the cluster 
(e.g., at their pericenter passages of orbital evolution), they 
can be identified as GGCs in the present selection method of 
ICGCs. Therefore we suggest that Egc of ICGCs in the cen- 
tral region of the cluster can be somewhat underestimated 
in Figure 2. 

The number fraction of ICGCs ranges from ~ 0.2 to 
~ 0.4 for the simulated clusters with 1.0 x 1O 14 M < M h < 
6.5 x 1O 14 M0. The power-law slopes of Egc range from ~ 

— 1.5 to ~ —2.5 for the clusters with the above mass range. 
About 10 — 20% of all GCs are located in the central 50 kpc 
of the simulated rich clusters, which implies that these inner 
GCs formed at high redshifts (z > 6) can be responsible for 
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Figure 3. Distributions of R ne i (left) and ikf h nei (right) for 
INGCs (solid) and for HGCs (dotted). Here R ne i represents the 
distance between a GC and the center of a halo that is nearest 
to the GC, and Mh jne i is the mass of the nearest neighbor halo. 
Therefore, R ne i and Mh >ne i for HGCs are those of their host halos. 




Figure 4. Distributions of R/Rh,gc (left) and (right) for 
INGCs (solid) and HGCs (dotted) at z = 6. Here R, R nc[ , and 
M n represent the distance of a GC from the center of its host 
halo at z = 6, the half-number radius of the GC system in the 
host halo, and the mass of the host halo, respectively. 



high Sn of the central cD galaxies. The more details on the 
parameter dependences of ICGC properties will be discussed 
in our forthcoming papers (BY). 



3.2 INGCs 

INGCs in the fiducial model comprises about 1% of all GCs 
formed in subhalos before z — 6 (See Figure 1). Although 
these INGCs are formed as a result of tidal stripping of 
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Figure 5. Dependences of the number fraction of INGCs 
(/ingc) on the truncation radii of GCSs at z = 6 (Rtr,gc) for 
the models with initial GC profiles similar to density profiles of 
dark matter halos (filled squares) and power-law ones with the 
slope of —3.5 (open squares). For comparison, the model in which 
only the very central particle in each subhalo is regarded as a GC 
(represented by "nuclei") is shown by a filled circle. 



GCs from subhalos during hierarchical structure formation 
through interaction and merging of subhalos between z — 
and 2 = 6, there appears to be no remarkable differences 
in the large scale distributions between INGCs and HGCs 
at 2 = (See Figure 1). Figure 3 however indicates clear 
differences in 7? ne i and M hl nci distributions between INGCs 
and HGCs. For example, the number fraction of INGCs with 
i?nei > 1 Mpc is 0.46 whereas that of HGCs with R ncl > 1 
Mpc is 0.09. The derived higher fraction in INGCs strongly 
suggests that INGCs are truly "free-floating" GCs in in- 
tergalactic/group/cluster regions. The Mh, nc i distribution 
of INGCs shows two peaks around M hl nci ~ 1O 11 M and 
« 1O 13 M and the number fraction of M h , nci > 1O 13 M 
is 0.27 for INGCs, which is significantly smaller than that 
(0.61) for HGCs. These results imply that the fraction of 
intercluster GCs among all INGCs can be observed to be 
small. 

Figure 4 describes (1) where progenitor GCs of INGCs 
were located with respect to the centers of their host sub- 
halos at 2 = 6 and (2) what the mass ranges of their host 
subhalos were at z = 6. The number fraction of GCs in the 
outer parts (R/Rh, S c > 2) of their host subhalos at z = 6 
is 0.17 for INGCs and 0.06 for HGCs in Figure 4. This re- 
sult suggests that GCs in the outer parts of subhalos, where 
GCs are more strongly influenced by external tidal force, are 
more likely to become INGCs and thus confirms that tidal 
stripping of GCs during interaction and merging of subhalos 
between 2 = and z = 6 is a major mechanism for INGC 
formation. The mean masses of host subhalos of INGCs and 
HGCs at 2 = 6 are 2.5 x 1O 1O M and 1.8 x lO n M , respec- 
tively, in Figure 4. This is also consistent with the above 
formation process of INGCs, because less massive subhalos 
are more strongly influenced by tidal stripping. If there are 
negative metallicity gradients of GC systems and positive 
relations between halo masses and GC mean metallicities in 
subhalos at z = 6, as observed for nearby GC systems (e.g., 
Ostrov et al. 1993; Cote et al. 1998), the results shown in 
Figure 4 imply that INGCs can be significantly more metal- 
poor than HGCs. 

Figure 5 describes how the number fraction of INGCs 
among all GCs (/ingc) depends on initial density profiles of 
GCSs at 2 = 6. It is clear from this figure that (1) /ingc is 
higher for larger Rti, gc and (2) this Rti, gc dependence can be 
seen both in the NFW profile and the power-law one with 
the slope of —3.5 (i.e., the observed profile of the Galactic 
GCS). Figure 5 also shows that only ~ 0.5% of GCs that 
were initially in the nuclear regions of subhalos at z = 6 can 
finally become INGCs. These results imply that /ingc can 
range from ~ 0.005 to ~ 0.02 for a reasonable set of model 
parameters on initial density profiles of GCSs. 

/ingc also depends on the methods to identify halos 
and GCs. /ingc is 0.0121, 0.0075, 0.0043, and 0.0024 for 
Nmi n =10, 20, 50, and 100, respectively. The number frac- 
tion of INGCs that were nuclei at z = 6 (/ingc.n) are 0.0046, 
0.0014, 0.0003, and 0.0001 for 7V min =10, 20, 50, and 100, re- 
spectively. This /ingc.n might well depend on the resolution 
of simulations, /ingc is 0.0067 in the model with the FoF 
linking length of 0.025. About 80 % of GCs that are identi- 
fied as INGCs are not bounded gravitationally by any closest 
halos. These dependences imply that /ingc can range from 
an order of 0.1% to 1%, given some uncertainties in the best 
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possible parameter values of the methods (e.g., FoF linking 
length) . 



4 DISCUSSIONS AND CONCLUSIONS 

Although van den Bergh (1958) already suggested the ex- 
istence of INGCs in the local universe, the present study 
first showed that (1) INGCs can be formed during hierar- 
chical formation processes of galaxies and clusters, (2) these 
INGCs are about 1 % of all GCs formed before z — 6, 
and (3) they can be typically more metal-poor than those 
within virialized galaxy-scale halos at z — 0. These INGCs 
are highly unlikely to suffer destruction processes by strong 
galactic tidal fields that are suggested to be important for 
understanding the origin of the observed mass function of 
GCs (e.g., Fall & Zhang 2001). Therefore, we suggest that 
(1) the mass function of INGCs can be significantly different 
from that of GCs in galaxies and (2) INGCs possibly retain 
fossil information on GC mass function at the epoch of GC 
formation. 

Then how many INGCs are expected to be observed in 
the intergalactic regions close to the Galaxy ? We can pro- 
vide an answer for this question by using the present result 
on the number fraction of INGCs and the initial GC num- 
ber of the Galactic GCS before GC destruction. McLaughlin 
(1999) showed that total number of initial GCs in a galaxy 
can decrease by a factor of 25 within the Hubble time ow- 
ing to GC destruction by the combination effect of galactic 
tidal fields and internal GC evolution (e.g., mass loss from 
massive and evolved stars). This means that the initial GC 
number is about 4000 for the Galaxy with the observed GC 
number of 160 at z = (van den Bergh 2000). By using 
the present result that ~ 1% of all GCs can become INGCs, 
the expected number of ICGCs in the intergalactic regions 
close to the Galaxy can be ~ 40. We thus suggest that some 
bright objects of these ~ 40 intergalactic GCs can be found 
in currently ongoing "all-objects" spectroscopic surveys for 
targeted areas (e.g., 6dF). 

Since White (1987) pointed out that clusters of galaxies 
might contain ICGCs, several authors have suggested some 
observational evidences for or against the existence of ICGCs 
(e.g., West et al. 1995; Blakeslee 1997; Harris et al. 1998). 
Our future more extensive numerical studies of ICGC forma- 
tion will provide observable predictions on some correlations 
between number of ICGCs and global properties (e.g., mass 
and X-ray temperature of hot gas) in their host clusters and 
thus help observers to confirm the existence of ICGCs. Fur- 
thermore, if physical properties of ICGCs strongly depend 
on Ztorm after which GC formation was severely suppressed 
by some physical processes (e.g., cosmic reionization) , they 
can provide some observational constraints on Zf orm . 
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